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Abstract. We summarize the main results obtained recently by our group on the identification and 
study of very high-z galaxies (z>7) using lensing clusters as natural gravitational telescopes. A descrip- 
tion of our pilot survey with ISAAC/VLT is presented, aimed at the spectroscopic confirmation of z>7 
candidate galaxies photometrically selected from deep near-IR, HST and optical ground-based imaging. 
The first results issued from this survey are discussed, in particular the global photometric properties of 
our high-z candidates, and the implications for the global star formation rate at very high-z. 



1 Introduction 

Considerable advances have been made during the last decade in the exploration of the early Universe, from the 
discovery and detailed studies of redshift z ^ galaxies (the so-called Lyman break galaxies, LBGs, e.g. Steidel 
et al. 2003), over z 4-5 galaxies found from numerous deep multi- wavelength surveys, to galaxies at z ~ 6-7, 
close to the end of reionisation epoch of the Universe (e.g. Hu et al. 2002, Kodaira et al. 2003, Cuby et al. 2003, 
Kneib et al. 2004, Stanway et al. 2004, Bouwens et al. 2004b). Extending the searches beyond z ~ 6.5 and 
back to ages where the Universe was being re-ionized (cf. Fan et al. 2002) requires extremely deep observations 
in the near-IR bands. Indeed, astounding depths can be reached in ultra-deep fields, such as demonstrated e.g. 
recently with J and H imaging of the NICMOS Ultra-Deep Field (UDF; Thompson et al. 2005; Bouwens et 
al. 2004a, Bouwens et al. 2005) from which 5 faint {Hab ~ 27) candidates at z ~ 7-8 have been identified 
(Bouwens et al. 2004b). 

We present in this paper a summary of our results on a deep survey of lensing clusters with ISAAC/VLT, 
aimed at constraining the abundance of star- forming galaxies at z '--^ 6 — 11 taking benefit from lensing magnifi- 
cation to improve the search efficiency and subsequent spectroscopic studies (see more details in Pello et al. 04 
and Richard et al. 05). We briefly describe the photometric technique used to identify very high-z objects, the 
construction and analysis of the photometric catalogs, the luminosity funtions (LFs) derived for star-forming 
galaxies up to z ~ 10, and the implications for the cosmic SFR. Throughout this paper we adopt standard 
cosmological parameters (JIa = 0.7, flm = 0.3, Hq = 70 km s~^ Mpc~^). 



2 Photometric Survey and Selection of High-z Candidates 

The objective of our Survey was to obtain deep near-IR photometry from 1.0 to 2.4 /im, in order to derive 
accurate photometric redshifts for optical dropouts in the critical domain 6 < z < 11. Simulations have 
been done to dcflne the observing strategy to target high-z sources using the evolutionary synthesis models 
by Schaerer (2002, 2003) for Population III and extremely metal deficient starbursts, together with the usual 
templates for normal galaxies. The main relevant signatures of genuine star-forming sources at z > 7, which are 
common to all models, are well known: they are optical dropouts, displaying a strong break and "red" optical 
vs. IR colors, whereas they exhibit a "blue" SED redwards due to UV rest-frame emission. Different redshift 
intervals are defined using the appropriate set of near-IR filters in combination with optical data through the 
Lyman break technique, by constraining both the position of the break and the restframe UV slope. Figure 1 
illustrates this technique in the 8 < z < 11 domain. 

The two lensing clusters used in this study were AC114 (z = 0.312) and Abell 1835 (z = 0.252). AC114 is 
a well-known gravitational telescope, with a lens model well-constrained by a large number of multiple-images 
at high-z (Small et al. 1995, Natarajan et al. 1998, Campusano et al. 2001). A1835 is the most X-ray luminous 
cluster in the XBACS sample (Ebeling et al. 1998), thus potentially one of the most efficient gravitational 
telescopes. We observed these clusters in the ~ 0.9 to 2.2 ^m domain between September 2002 and April 2004, 
covering as far as possible the z, SZ, J, H, and K bands. Optical images between U and / bands were available 
from previous surveys and data archives (Table 1). 

Near infrared photometry of extremely faint sources requires a careful data reduction, described in details 
by Richard et al. (2005). In summary, after a standard pre-reduction (ghost, dark and flat-field corrections). 



Table 1: Summary of the photometric dataset used in this Survey: filter, exposure time, seeing (original images), 
pixel size, 1 tr limiting magnitude within 1.5 " diameter aperture, filter effective wavelength, AB correction 
(toab = rnvega + C ab) ^ and references (see Richard et al. 2005 for details). 
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we used the IRAF package XDIMSUM^ for a two-step sky-subtraction. During the first pass, each image is 
sky-subtracted using the sky pattern obtained from a group of adjacent frames and a bad-pixel mask is created 
in the process. Images are registered and combined using integer shifts values to preserve the noise properties, 
with bad-pixel rejection. Then, sources are detected in order to create an object mask, and a second sky- 
subtraction is applied to the data. Several versions of the final images were produced, using slightly different 
reduction recipes, in order to cross-check the final catalogs. 

Photometry in the near-IR bands was obtained after matching all images to a common seeing with a gaussian 
convolution, the worst case being the J band for both clusters. The SExtractor package (Berlin & Arnouts 1996) 
was used for source detection and photometry. We optimized the parameters to detect very faint unresolved 
sources, in order to build an i?-band selected sample. The original images were used to derive the error bars in 
each band through detailed simulations. Also limiting magnitudes in the Table 1 and completeness levels in the 
different filters were obtained in this way. The final catalogs include photometry within 1.5" apertures for all 
objects detected in the H band; we were able to measure photometry of very faint sources {J ~ 24.4 — 24.8, H 
and Ks ~ 23.5) with a relatively good accuracy (S/N> 3-4). The fraction of spurious detections expected in our 
photometric catalogs, for objects detected only in the reference filter H, was estimated from a special H band 
noise image where all astronomical sources were removed by subtracting by pairs sequential images acquired 
with similar seeing conditions, and then coadding them using the same procedure as for astronomical images. 
The result is an image with the same noise properties compared to the final stacks, but without astronomical 
sources. We find that for the faintest magnitude bins considered in this survey, the maximum number of false- 
positive detections should be typically lower than 50 %, for sources detected only in the reference filter H, and 
no spurious detections are expected up to = 23.0 (AB~ 24.4). 

A catalog of optical dropouts was selected (i.e., objects non-detected in all the optical images, from U to z 
bands). First-category sources are those detected in at least two near-IR bands. Among them, ~ 89% of sources 
detected in the iJ-band reference image in A1835 (and 75% in AC114) are also re-detected in the pseudo-x^ 
image. The later was obtained from the combination of individual J, H and Ks images, normalized by the 
noise la image, and weighted by the square root of the corresponding exposure-time maps. After careful manual 
inspection, the final catalog of first/second-category dropouts for Abell 1835 and AC114 contains 18/6 and 8/2 
sources respectively, and number of third category sources (dubious after manual inspection, detected only in 
the reference filter) close to the critical lines. 

The position of optical dropouts in the different color-color diagrams provides an estimate of their photo- 
metric redshift, and an objective criterium to classify them into different z intervals. The majority of optical 
dropouts in the two clusters fulfill the high-z requirements in Fig. 1. Those located in the < z < 8 region 
of the diagram fulfill the ERO selection criterium {R — K > 5.6), and some of them could be intermediate- 
redshift dusty starbursts. The SZJH and zSZJ color-color diagrams were used to select candidates in the 
range 7 < z < 8.5 and 6 < z < 7.5 respectively. For about 30% of our candidates, the S/N is enough to derive 
photometric redshifts based on SED-fitting using an adapted version of the public software Hyper z (Bolzonella 
et al. 2000). Objects unambiguously identified as low-z galaxiess are excluded from the sample, as well as 
sources whose nature could not be determined with the present data (either variable sources or "bright" EROs 
with ambiguous but plausible low-z solutions). 

^XDIMSUM is a modified version by the IRAF group of the Deep Infrared Mosaicing Software package by P. Eisenhardt et al. 
See ftp://iraf.noao.edu/extern-v212/xdimsum for details 
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Figure 1: a: J — H versus H — Ks color-color diagram (Vega system) showing the position expected for different 
objects over the interval z ~ to 11. The position of stars and normal galaxies up to z < 8 are shown, as 
well as the shift direction induced by Ay = 1 magnitude extinction. Thin and thick lines display models below 
and above z = 8 respectively. Several models for Pop III starbursts are presented, for different fractions of 
Lyman-a emission flux entering the integration aperture: 100% (red solid line), 50% (red dashed line) and 0% 
(red dot-dashed line). The location of Kinney et al. (1996) starbursts templates is also given for comparison 
(SBl(cyan) and SB2 (blue)). All star-forming models enter the high-z candidate region at z > 8. b: Location of 
first and second category photometric candidates in the lensing cluster AC114 (circles and crosses respectively). 
Contours are overplotted for magnification values of 1, 2 and 3 magnitudes, computed assuming sources at 
z = 9, although the position of these lines is weakly sensitive to source redshift within the relevant z 6 — 10. 

High-z candidates were selected based only on their photometric properties, irrespective of their positions 
with respect to the critical lines. However, objects located close to the high-z critical lines are of greater interest, 
because of the larger magnification. The minimum magnification factor over the region covered by our near-IR 
survey is at least ~ 0.7 magnitudes, and at least ~ 1 magnitude over 50 % of the ISAAC field of view. Thus, 
the effective 3ct limiting magnitudes reached here are, at worst, similar to those attained in the HDFS (Labbe 
et al. 2003) in JHKs (respectively AB^^ 26.8, 26.2 and 26.2). Our 3a limiting magnitudes in the H band are 
also very close to the typical magnitudes of the z ~ 7 — 8 z-dropouts detected by Bouwens et al. (2004b) in the 
Hubble Ultra Deep Field, with Hiqq{AB) ~ 26.0 to 27.3, after correction for a typical magnification factor of 
at least ^ 1 magnitude. 



3 Luminosity Function and Cosmic Star Formation Rate 

The typical magnification values of our candidates range between 1.5 (^ 0.44 mags) and 10 (2.5 mags). For 
some objects very close to the critical lines, we found magnifications values /i > 25. Interestingly, although 
the selection criteria are only based on near-IR colors irrespective of magnitudes, almost all the photometric 
candidates fulfilling our selection criteria turn out to be fainter than H = 23.0 (AB ~ 24.5). Only three 
exceptions are found in Abell 1835 among the possible low-z EROs, as described above. After correction for 
magnification accross these fields, the lack of "bright" sources means that we have not detected young starbursts 
at z ~ 6 — 10 more massive than typically a few 10* M0 (under standard assumptions for the IMF). 

We derived the unlensed ii5oo luminosity, at 1500 A restframe, for all high-z candidates, using the adopted 
photometric redshift. L1500 luminosities were converted into Star Formation Rate (SFR) through the usual 
calibration from Kcnnicut (1998). The typical SFR obtained for objects included in the final sample is ^ 10 
Mq yr^^, with extreme values ranging from a few units to ~ 20 Mq yr~^. The restframe UV slope of our 
candidates is extremely blue, usually ranging between —1.5 and —3.5, a systematic trend also reported by 
Bouwens et al. (2004b) for their sample of z ~ 7 — 8 candidates. Although optical-dropouts are stretched by 
the magnification factor they appear as point-like sources in our ground-based images. The physical size of 
these objects at z > 7 is likely to be smaller than 1.7 kpc, with the magnification factors involved. 

Magnification and dilution effects by the lensing field were carefully taken into account to compute number 
densities and derived quantities, in particular to estimate the LF at 1500 A . The observed sample of candidates 
was also corrected for incompleteness using mock simulations (see details in Richard et al. 2005). The combined 
L1500 LFs for both clusters, with the corresponding error-bars, are given in Fig. 2. Only first-priority candidates 




Figure 2: a: ^1500 LFs derived for the photometric sample of high-z candidates in this survey (adapted from 
Richard et al. 05). z intervals shown are: z = 6 — 10 (red) and z = 8 — 10 (black). Errors include Poisson noise 
statistics and lensing uncertainties. STY fits to the LP data are presented by thick solid lines, and compared 
to the LP by Steidel et al. for LBGs at 2 4 (thick dashed line), with the usual (1 + z)~^ correction, without 
any additional renormalization. We also display the LFs corresponding to simple models, for 2 extreme IMP 
assumptions discussed in the text, b: Evolution of the comoving SPR density as a function of redshift. Our 
upper limits are compared to other surveys, uncorrected for extinction (adapted from Bunker et al. 2004, and 
references herein): data compiled from the GPRS (Lilly et al. 1996, filled triangles), Connolly et al. 1997 (filled 
pentagons), LEG work from Steidel et al. 1999 (open squares), Pontana et al. 2003 (open circles), Iwata et 
al. 2003 (cross), Bouwens et al. 2003a (filled diamonds), GOODS (Giavalisco et al. 2004, open stars), ACS 
estimates by Bouwens et al. 2003b (filled stars) and Bunker et al. 2004 (filled circle). The value derived by 
Bouwens et al. (2004b) in the UDP is also shown. Our results are displayed for all (solid) and first-category 
(dashed) candidates. 



have been considered, but the difference obtained when using the full sample is within la error bars. STY fits 
(Sandage, Tammann & Yahil 1979) to the data are also presented in Fig. 2. The typical value found for L* 
is 10^^'^ ergs s~^ s~^ A with a fixed value a = 1.6 (i.e., Steidel et al. (1999) determination for LBGs at 
z ~ 4). The STY fit to the data is in remarkably good agreement with the LP found by Steidel et al. for LBGs 
at z ~ 4, with the usual (1 + z)~^ correction to account for the surface brightness increase with redshift due 
to size scaling for a fixed luminosity, without any additional renormalization. A fairly good agreement is also 
found when comparing with the LP derived by Bunker et al. (2004) for their sample of z ~ 6 candidates in the 
UDP, i.e. a density of sources ~ 1/6 (-0.8 dex) smaller as compared to LBGs at z ~ 4, still within our Ict error 
bars. We compare the observed LFs to the predictions obtained from a simple model for halo formation based 
on Press-Schechter formalism, assuming that all haloes convert a constant fraction ~ 0.1 of their baryonic mass 
into stars somewhen between z = 17 and z = 6, with a correction for the visibility time of starbursts, and using 
two extreme IMP assumptions: a standard Salpeter and a top-heavy IMP (stars between 50 and 500 Mq). "Top 
heavy" IMF models provide a better fit for the bright end of the LP, but this simple model can hardly explain 
simultaneously the behaviour of the bright and faint ends of the LP. 

Figure2 displays the upper limits for the Cosmic SPR value obtained for each redshift bin, by integrating 
the LFs down to 0.3 ^^=3, compared to other surveys. 



4 Discussion and Conclusions 

Taken at face value, the cosmic SPR density found in this survey is in good agreement with the theoretical 
estimates for the redshift domain considered here derived by Barkana & Loeb (2001; see their Pig. 29), for a 
reionization redshift ~ 8 — 6. However, there is a discrepancy by a factor of ~ 10 between our results and previous 
studies at similar redshifts, in particular in the UDP (Bouwens et al. 2004, 2005). In all cases, the sources 
detected are photometric candidates, and thus upper limits to the actual UV flux densities. The effective fields 
surveyed are dramatically small, thus leading to strong field-to-field variations in the number of sources. Gluster- 
to-cluster fiuctuations are clearly seen in our sample, although lensing and photometric considerations could 
account for most of them. A positive magnification bias could still be present in this survey (the incompleteness 
of our sample in the relevant magnitude domain is smaller than in blank-field surveys), producing a systematic 
trend as compared to blank fields. Could this result be confirmed on a larger sample of lensing clusters and 



blank fields, the slope of the number counts at z ^ 6 — 10 could be precisely constrained, at least for the 
brightest part of the LF. 

Up to now, our spectroscopic survey with ISAAC has targeted 2 candidates in AC114, and 7 in Abell 1835 (4 
first priority targets and 3 secondary ones; Pello et al. 2004); 2/3 of objects in this sample display emission lines. 
The efficiency of our survey nowadays could range between 30 and 50%, with interesting low-z by-products. 
A large majority of our high-z candidates still need to be confirmed, either by a redetection of the faint emission 
line, or by the non-detection of other lines expected at low-z. 

The results presented here are to be confirmed in different ways. An enlarged spectroscopic survey is urgently 
needed to determine the efficiency of our selection technique. Also, increasing the number of lensing fields with 
ultra-deep near-IR photometry is essential to get tighter constraints on the abundance and physical properties 
of z > 7 starburst galaxies. 
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